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0\ • ABSTRACT 

The claim by Dar and Shaviv that they have found a standard model solution 
to the solar neutrino problem is based upon an incorrect assumption made in 
extrapolating nuclear cross sections and the selective use of a small fraction of 
the nuclear physics and of the neutrino data. In addition, five different solar 
model codes show that the rate obtained for the chlorine experiment using the 
Dar-Shaviv stated parameters differs by at least lAa from the observed rate. 
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1. Introduction 



In a widely circulated preprint, Dar and Shavivcl claim to have provided a 
standard solar model solution to the solar neutrino problem. Their preprint has 
aroused sufficient curiosity that we feel compelled to point out its most obvious 
defects, although we do not plan to publish our response in a refereed journal 
unless the Dar and Shaviv preprint is also published in a refereed journal. 

To summarize our response, Dar and Shaviv have not solved the solar neutrino 
problem. They have made an incorrect assumption in extrapolating nuclear cross 
sections and have used nuclear and neutrino data selectively. Most surprisingly 
their solar model results are not reproducible. We obtain from five independent 
solar neutrino codes and the Dar-Shaviv input parameters, neutrino fluxes that 
are inconsistent with the values reported by Dar and Shaviv. With the Dar 
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and Shaviv input parameters, all four solar neutrino experiments differ from the 
solar model results (obtained with the five different codes) by more than 14a 



for the chlorine experiment 



more than 



^Rate (theory) — Rate (exp)^ /cr(exp) 

4a for the Kamiokande experiment,! more than 4a for the GALLEX (gallium) 
experiment,Q and more than 3a for the SAGE (gallium) experiments 

This response is organized as follows. In §1, we comment on the fact that Dar 
and Shaviv chose to compare their results with only a small part of the available 
chlorine solar-neutrino data. In §2, we point out that they made an ad hoc and 
incorrect assumption in extrapolating the nuclear physics data and compared 
their results with only a small fraction of the available experiments. In §3, we 
show that their solar model calculations are in disagreement with calculations 
made using five other well-tested codes. 

2. Choosing Part of the Chlorine Solar-Neutrino Experimental 
Data 

Dar and Shaviv chose to consider (see their Figure 1) only four years of data 
from the chlorine solar neutrino experiment beginning in 1987, although 23 years 
of data have been reported 2 . For the period beginning in 1987, the experimental 
measurement of the chlorine rate is 2.8 ± 0.3 SNU. The measured rate for the 
entire period for which data has been reported is 2.28 ±0.23 SNU. The measured 
rate in the chlorine experiment, during the short period considered by Dar and 
Shaviv, is still more than 4a less than their calculated result of 4.2 SNU and is 
more than 14<t from the value we calculate (see §3 and Table 1 below) for their 
stated parameters. 

3. Extrapolating Nuclear Cross Sections with an ad hoc As- 
sumption and Using Only Part of the Experimental Data 

Dar and Shaviv claim that by replacing the conventional! representation of 
the extra energy dependence (in addition to the point-nuclei barrier penetration 
factor) of the nuclear cross section factor by another smoothly-varying, param- 
eterized function, they obtain a very different extrapolation to low energies for 
the rate of the 3 He(a,^/) 7 Be reaction. Of course, both procedures must yield the 
same answer if they are both applied in an unbiased way to the same experi- 
mental data set and all the energy dependences are included (see below). For the 
reactions discussed by Dar and Shaviv, the existing data sets determine well the 
extrapolations. 

Dar and Shaviv, however, assumed incorrectly that the only energy-dependent 
effect besides point-nuclei barrier penetration is nuclear size, and then compared 
their results with only two of the nine published experiments. Their selective use 
of data and this incorrect assumption explain why the Dar and Shaviv answer for 
S(0) differs from the standard value obtained by nuclear physicists. We explain 
in more detail below. 

For extrapolation of nuclear cross section data to stellar energies, it is con- 
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ventional to define the astrophysical S-factor, 



S(E) = a{E) E exp{ 2nr)(E)}. (1) 
The Sommerfeld parameter is given by 

where v is the relative velocity of the two particles in the entrance channel, 
E is their relative energy, and Z\,Z<i are their charges. The form of Eq. (1) 
accounts explicitly for the energy dependence of s-wave tunneling through the 
Coulomb barrier of two pointlike particles and a kinematic flux factor. In the 
absence of near-threshold resonances, the energy dependence of the S-factor is 
expected to be small at low energies. However, in order to carry out a reliable 
extrapolation, the energy-dependent effects that are not accounted for in Eq. (1) 
must be allowed to show up in S(E). These effects include nuclear structure, the 
strong interaction, energy dependent operators in the transition matrix elements, 
antisymmetrization between the colliding nucleons, finite nuclear size, the final- 
state phase space, and the contributions from other partial waves. In fact, for 
all the reactions of importance for the solar p — p chain, the observed energy 
dependence of the various reactions agrees well with that calculated in theoretical 
models which account explicitly for the known nuclear effects that are omitted 
from Eq. (1). In particular, it was demonstrated more than thirty years ag 
that the complete Coulomb wave function provides a good description of the 
measured cross section factor for the 2 H(p,^/)~ i He reaction down almost to solar 
thermal energies (16 keV). 

Dar and Shaviv have chosen to factor out a slightly different energy depen- 
dence from the cross section data by defining a modified S-factor, S(E), as 

S(E) = a(E)Eexp{2irr](E)[l + 2/vr(2 v / ^ - arcsin^/x - yjx(l - x)))}, (3) 

with x = E/E c and the Coulomb energy E c = Zi^e 2 /R at a radius parameter 
R. With this definition, S(E) attempts to account for the finite size of the nuclei. 
However, S(E) is still expected to be dependent on E, because of the other effects 
listed above that introduce an energy dependence in addition to the finite size 
of the nucleus. The theoretical models that have been used previouslyH^ to 
extrapolate the cross section data to solar energies properly take account of the 
finite nuclear size effects along with the other effects discussed above. Obviously 
both definitions, Eqs. (1) and (3), must lead to the same results for the low energy 
S-factor when proper account is taken of the additional energy dependence not 
included explicitly in the respective equations. 

Dar and Shaviv did not take account of the additional energy dependences 
nor of all of the available nuclear physics data. For the determination of the low- 
energy cross section factor, S34, for the 3 He(a, 7) 7 Be reaction, Dar and Shaviv 
apparently adjusted the radius parameter R (see Eq. 3 above) so that the energy 
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dependence of S(E) is mostly removed for two of the nin existing experiments. 
(They seem not to have noticed that the value of R that they obtain is very 
different from the measured radius of 2.8 fm determined by electron scattering.) 
They did not allow for the other energy dependences discussed above and they 
only took account of two of the experiments. 

For the determination of the cross-section factor for the 7 Be(p, 7) 8 -B reac- 
tion, S17, Dar and Shaviv first remove the contribution from d-wave capture 
(based on previous calculations of this single component of the energy depen- 
dence). However, the d-wave capture contributes about 6% to the S-factor at 
the most effective energy for the reaction (about 20 keV) and nmst be added to 
the computed s-capture rate. In addition, previous calculationaj'EJ demonstrate 
that there are significant further contributions to the energy dependence for S17 
that must be included in order to make a reliable extrapolation. When these 
additional effects are included in the extrapolation of Sn(E) (Eq. 3 above), the 
results must agree with the standard extrapolations B 

Dar and Shaviv cited the preliminary Coulomb dissociation work described 
in preprint formll3 as evidence for a lower-than-standard value for the crucial 
cross section factor for the 7 Be(pjx) s B reaction. When the E2 contribution to 
this reaction is taken into account ,£3 the preliminary Coulomb-dissociation value 
differs from the six direct measurements^ of the 7 Be(p,~/) S B cross section by 
a factor of two while the estimated uncertainty^ in direct measurements is only 
11%. Moreover, there are still some unanswered questions about the application 
of the Coulomb-dissociation method for determining radiative capture cross sec- 
tions, aside from the experimental difficulties inherent in covering a sufficient 
range in energy and angle to validate the reliability of any inferences. 

For other nuclear reactions, Dar and Shaviv have used the low-energy cross 
section factors from an earlier revie w0 which provided fitting formulae suit- 
able for use at temperatures (~ 10 9 K) much higher than are reached in the sun 
(~ 10 7 K). The quantitative effect of these approximations is difficult to estimate, 
especially since other authors (see section 4) use-for solar calculations-explicit 
formulae that are suitable for the lower solar temperatures. However, an approx- 
imate discussion of using the fitting formulae for higher temperatures at solar 
temperatures has been givenll3; the principal effects of the high-temperature for- 
mulae are in the direction to decrease the predicted 7 Be and 8 B neutrino fluxes. 

Dar and Shaviv discuss at some length the fact that the Debye-Hiickel ap- 
proximation to the screened nuclear potential is not correct everywhere in the 
sun. It is not clear what they recommend (although they say the effect is small), 
nor if they are aware that modern screening calculations go well beyond what 
they discuss. E30 

4. Solar Model Calculations 

Six authors of this paper (Bahcall, Christensen-Dalsgaard, Degl'Innocenti, 
Glasner, Pinsonneault, and Proffitt) have repeated the solar model calculations 
of Dar and Shaviv using the non-standard parameters that Dar and Shaviv chose, 
namely, a solar luminosity of 3.826 x 10 33 erg s _1 and low energy cross-section 
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factors of 534(0) = 0.45 kev-b and £17(0) = 17 eV-b. Dar and Shaviv did not 
specify in their preprint many of the important input quantities in their model; 
they did not state what they used for the element abundances, the radiative 
opacities, the equation of state, and the neutrino cross sections. They did not 
say which of the several available perscriptions for diffusion they used. We have 
therefore carried out calculations using a variety of different choices for these 
quantities, namely, the choices made previously as their best estimates by the 
six different authors who used five independent stellar evolution codes .EJiH3 

The results are shown in Table 1. The first column identifies the computer 
code used in constructing the solar model; the second column indicates whether 
or not particle diffusion was included. The third, fourth and fifth columns give the 
calculated rates predicted for the chlorine, gallium, and Kamiokande solar neu- 
trino experiments, respectively. The numbers in parentheses indicate the number 
of standard deviations quoted by the experimentalists (adding statistical and 
systematic errors quadratically) by which the calculated rates differ from the 
measured rates. For gallium, we have compared with the most recent GALLEX 
determination (which has smaller quoted uncertainties). The last two columns 
of Table 1 give the fluxes of 13 N and 15 solar neutrinos. Some details regarding 
the solar models are given in the footnotes to the tabic. 

The calculations reported by Dar and Shaviv, the first row of Table 1, predict 
a rate that differs by 8a from the average chlorine experimental value. Even if 
their calculations and assumptions were correct (and we think they are not), then 
they have not solved the solar neutrino problem unless one completely disregards 
the chlorine solar neutrino experiment. 

In rows two through four we give results obtained with five different computer 
codes that have been intercompared in the literature with other standard models; 
all of these codes have previously been shown to give the same results for the 
same input parameters and to give results that agree well with still other well- 
tested solar model codes. For the Dar and Shaviv parameters, all of the five 
well-calibrated solar codes give results very close to each other, as they must, 
but they all differ significantly from the Dar and Shaviv results. The relatively 
small differences between the five well-calibrated solar models are primarily due 
to the treatment of particle diffusion (see column two). This is a significant 
physical effect that is somewhat difficult to include in solar model calculations, 
but which has been taken into account in the most recent and advanced solar 
model computer programs. 

For the Dar and Shaviv parameters, the five well-calibrated models give for 
the chlorine experiment between 5.6 SNU (14a discrepancy) and 6.3 SNU (17a 
discrepancy) depending upon whether or not particle diffusion is included. This 
contrasts with the value of 4.2 SNU reported by Dar and Shaviv. All five of the 
models predict about 125 SNU (4a discrepancy) for gallium, which contrasts 
with the Dar and Shaviv value of 109 SNU. Finally, the 8 B flux in the five 
well-calibrated solar models is between 4a and 6a from the Kamiokande result 
(depending upon whether or not particle diffusion is included), while the Dar 
and Shaviv result agrees almost exactly with the Kamiokande measurement. 
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Even more striking differences exist for the N and O neutrino fluxes, for 
which the Dar and Shaviv values differ by an order of magnitude from the values 
obtained in the other programs. Their flux for the 13 N neutrinos is low by about 
a factor of 6 and their 15 flux is low by about a factor of 15. 

None of the well-tested solar codes that we have used are able to reproduce 
the Dar and Shaviv results for neutrino fluxes. 

5. Summary 

Dar and Shaviv did not succeed in solving the solar neutrino problem despite 
the introduction of an ad hoc assumption, the selective treatment of data, and 
the incorrect extrapolation of some nuclear cross sections. Their failure to solve 
the problem is not surprising since it has been demonstrated elsewhereE2l that 
the chlorine and the Kamiokande experiments are inconsistent with one another 
if one assumes (on the basis of standard electroweak theory) that the shape of 
the 8 B neutrino spectrum is the same in the laboratory and in the sun. 

Several of the authors (JNB, SD, AG, KL, PDP, MHP, and CRP) thank the 
Institute for Nuclear Theory at the University of Washington for its hospitality 
and the Department of Energy for partial support during the completion of this 
work. 
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Table 1 

Solar Model Results for Dar and Shaviv Parameters 



Solar Code 


Diffusion 


37 C1 
(SNU) 


71 Ga 
(SNU) 


8 B t 


13 N * 


15Q* 


Dar and Shaviv 


Vac? 

i es i 


A 9^8 rr\ 




Z.o^U u J 


n 7 

U. 1 


n 9 
u.z 




ies 


C 9/17 rr\ 


izo^4 a ) 




K Q 

o.y 


0.1 


AARHUS b 


No 


5.6(14o-) 


122(4o-) 


3.6(4 a) 


5.4 


5.4 




Yes 


6.1(16o-) 


124(4 a) 


4.0(5 a) 


5.8 


5.8 


Promtt c 


Yes 


6.3(17 a) 


126(4 a) 


4.2(6 a) 


6.3 


5.5 


FRANEC d 


No 


5.8(15o-) 


123(4 a) 


3.9(5o-) 


4.9 


4.2 


ASTRA e 


No 


5.7(15o-) 


123(4 cr) 


3.7(4 a) 


4.2 


3.5 


T Units: ID 6 cm" 


V 1 ; *Units 


: 10 s cm _2 s 


-l 









a cf. Ref.E±l Includes helium and heavy element diffusion, as well as other improvements. 
b cf. Ref.0 Second model includes helium diffusion. 
c cf. RefH Includes helium and heavy element diffusion. 
d cf. Ref.0 Does not include diffusion. 

e cf. RefJnl Analagous to the Best model without diffusion of refill 
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